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Abstract 

We present analytical results for the biased diffusion of particles moving under a constant force in 
a randomly layered medium. The influence of this medium on the particle dynamics is modeled by 
a piecewise constant random force. The long-time behavior of the particle position is studied in the 
frame of a continuous-time random walk on a semi-infinite one-dimensional lattice. We formulate 
the conditions for anomalous diffusion, derive the diffusion laws and analyze their dependence on 
the particle mass and the distribution of the random force. 
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I. INTRODUCTION 



A vast variety of physical, chemical, biological and other natural processes can be ad- 
equately described by random processes exhibiting anomalous diffusion behavior at long 
times. This behavior, which is characterized by a nonlinear dependence of the variance of 
these processes on time, can be observed in various systems. Anomalous diffusion actually 
occurs, for example, in turbulent fluids amorphous solids [2], rotating flows [3], single 
molecules 41 and porous substrates jsl, and has been predicted to occur in many other 
systems 

The existence of anomalous diffusion in systems with quenched, i.e., time- independent 
disorder has also been extensively studied One of the most effective and simple ways to 

describe anomalous diffusion in these systems is based on the motion equations for diffusing 
objects (which we will call particles). In these equations, the influence of quenched disorder 
is usually modeled by a time-independent random potential producing the corresponding 
random force, and the influence of thermal fluctuations is accounted for by white noise. 
Specifically, this Langevin-type approach has been successfully applied to study a variety of 
phenomena, including biased diffusion, which occur when particles move under a constant 
external force in a one-dimensional potential 

If thermal fluctuations are absent then particles can be transported to an arbitrary large 
distance only if the distribution of the random force has bounded support. In this case, the 
directional transport of particles can be caused by either a periodic external force 12| or a 
constant one. In the latter case, particles move only in one direction and so the completely 
anisotropic case of biased diffusion, when the probability of motion along and against the 
external force equals 1 and 0, respectively, may exist. It has been shown for particular cases 
of the random force distribution that in the overdamped limit this diffusion is normal if the 
total force acting on a particle is strictly positive or strictly negative QQ. In Ref. Q 
it was also argued that the anomalous regimes of biased diffusion would exist if the lower 
(upper) bound of the total force at a fixed external force is equal to zero. However, none of 
the laws of anomalous diffusion was found in this case. 

The aim of this paper is to study the anomalous regimes of biased diffusion of particles 
moving under a constant force in a randomly layered medium which acts as a piecewise 
constant random force. The paper is organized as follows. In Sec. [Tll we describe the model. 
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reduce it to a continuous-time random walk (CTRW) on a semi-infinite chain, and calculate 
the first two moments of the particle position. The connection between the waiting time 
probability density and the particle mass and the probability density of the random force is 
also presented in this section. The conditions providing the anomalous behavior of biased 
diffusion are formulated in Sec. Illli In Sec. \IV\ using the Tauberian theorem and its modified 
version, we derive the laws of anomalous diffusion and analyze the influence of the particle 
mass and the random force distribution. Finally, in Sec. |V]we summarize our results. 

II. MODEL AND BASIC EQUATIONS 

We consider the one- dimensional propagation of a particle in a medium composed by the 
layers of a fixed width / whose transport properties are assumed to be random. The motion 
of a particle in this medium occurs under the action of a constant external force / (> 0), 
and the influence of the layers on the particle dynamics is modeled by a random force 
g{x). We assume that g{x) (i) is a bounded function, i.e., g{x) G [~go,go], (h) possesses a 
symmetry property, i.e., g{x) and —g{x) are statistically equivalent, and (iii) has statistically 
independent values on different intervals of the length /. In accordance with these conditions, 
we approximate g{x) by a piecewise constant random force (see Fig. [1]) whose values are 
distributed with the same probability density u{g). In this stage, we consider u{g) as an 
arbitrary symmetric probability density, u{g) = u{—g), satisfying only the normalization 
condition Jl°g^dgu{g) = 1. 

Since the total force acting on a particle equals / + g{Xt), where Xf (Xq = Xq = 0) is 
the particle position, its dynamics can be described by the motion equation 

^IXt + lyXt = f + g{Xt) (2.1) 

with II and u being the particle mass and the damping coefficient, respectively. According 
to this equation, ii f > go then / + g{Xt) > and particles can be transported to an 
arbitrary large distance in the positive direction of the axis x. Thus, in this case the 
condition Xt\^^^ oo holds for all sample paths of g{x). On the contrary, if f < go 
then for each sample path of g{x) there always exists a certain point L = Inst (^st = 
n = 0, 1, . . .), which is characterized by the conditions / > —g{x)\x<L and / < —g{x)\x=L+o, 
where particles are stopped. The probability Wn that Ugt = n is expressed through the 
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FIG. 1: Sample path of a piecewise constant random force g{x). 

probability / = jlg^ dgu{g) that g{x) < / as follows: Wn = 1^—1^+^, Therefore, the average 
distance (L) (the angular brackets denote an average over the sample paths of g{x)) from 
the origin to the stopping point can be written in the form (L) = I X^^i = I Z^^i 
Finally, using the geometric series formula, we obtain the desired result 

{L)=lj^. (2.2) 



Assuming that the probability that \g{x)\ = go equals zero, i.e., the probability density 
u{g) is not concentrated at the edges of the interval [—go, go], one can make sure that 1^1 
and so (L) — )■ oo as / — )■ (^o- Hence, at f = go the condition Xt\^^^ — )■ oo holds almost 
surely, i.e., with probability one. In contrast, if the probability density u{g) has unbounded 
support with go = oo as, e.g., for a Gaussian distribution, then J < 1 and so (L) is finite 
for all finite values of the driving force /. In other words, in this case particles cannot be 
transported to an arbitrary large distance. It is therefore we consider here only a class of 
probability densities u{g) with bounded support. It should be noted in this context that, 
since infinite values of g{x) are physically not relevant, the assumption of bounded support 
is not too restrictive. 

Our aim is to study the long-time behavior of the particle position Xt at f > go- The 
main statistical characteristic of Xt is its probability density function P{x,t) defined as 
P{x,t) = {6{x — Xt)), where 6{x) is the Dirac 6 function. If the solution of Eq. (12. ip were 
known for all sample paths of g{x), it would be, in principle, possible to determine P{x,t) 
directly from the definition. However, this approach is difficult to implement and, what is 
more important, it is not necessary for finding the long-time behavior of the moments of 
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Xf. Moreover, since at long times Xt can be accurately evaluated as a total length of the 
intervals {nl, nl + l) {n = 0,1, . . .) which a particle passes, many of the details of the particle 
dynamics described by Eq. fl2.ll) are needless for this purpose. 

It is therefore reasonable to consider, instead of the model based on Eq. fl2.ip . the uni- 
directional CTRW of a particle on a semi- infinite one- dimensional lattice with the period /. 
Introduced more than four decades ago [isl, the CTRW model has become one of the most 



effective and powerful tools in the theory of anomalous diffusion (see, e.g., Refs. |8NlO|). 
Within this model, we describe the particle position by a discrete variable Yt = lN(t), where 
N{t) is the random number of jumps up to time t. In order to guarantee that the long-time 
behavior of Yt and Xt are the same, we assume that for all n the waiting time r("), i.e., the 
time of occupation of the site nl, is equal to the time that a particle spends moving from 
the site nl to the site nl + l. If the inertial effects can be neglected then Eq. (12. ip yields 
r(") = z///(/ + 5f("^), where g^'^'^ = g(x) and x belongs to the nth interval, i.e., x E [nl,nl + l). 
Since the random forces Qn are statistically independent and distributed with the same prob- 
ability density u{g), the waiting times r*^"-* are also statistically independent variables whose 
probability density is given by 

— (— _ f\ \ 1 

p{r) = V ) , (2.3) 

0, otherwise 



where 



(2.4) 



J +90 f - 90 

It is important to emphasize that the inertial effects, at least in the underdamped regime 

characterized by the condition ut^i^/ ji ^ 1 (weakly underdamped regime), can also be 

incorporated into the CTRW framework. In order to illustrate this, let us first write the 

particle velocity v^'^\/S.t) [At G (0,r*^"'')] on the nth (n > 1) interval. The straightforward 

integration of Eq. ( 12. ip yields 

v^-\M) = ^±f^ + - l±f^y-^\ (2.5) 

where k = u/jj, and = v^^\0) is the particle velocity at the left end of the nth interval. 
Introducing also the particle velocity at the right end of this interval, ), from 

Eq. (12. 5p we obtain 



Then, taking into account that / = /J"'"^ d{At) v^^'^At), with the help of Eqs. (jZS]) and (Ej]) 
we find 

l=L±9_^(n)_Y± (2.7) 

U K 

Since in the case under consideration kt*^"^ ^ 1, the exponential term in Eq. (12 .Op can 
be neglected yielding = (/ + 9^^^)/^- According to this approximation, the particle 
velocity at the end of the nth interval is determined by the random force on this interval. 
Therefore, using the continuity condition for the particle velocity, f!?^ = ^\ we obtain 
^(") — (J _j_ Substituting these expressions for and f!"^ into Eq. (12.71) . we 

arrive to the following result: 

" «:(/ + ^?(")) ■ ^ ^ 

It shows that in the weakly underdamped regime the waiting time r^") depends not only 
on the random force g'^'^\ as in the overdamped case, but also on the random force g^"'~^\ 
Since these forces are statistically independent, the probability density of the waiting time 
can be written in the form 

p{r) = r r dgdg'u{g)u{g') 6 (r - ^^TT^l ^2.9) 

if r G [Tinin, Tinax]) otherwisc it equals zero. It is not difficult to verify that in the overdamped 
case (when k = 00) Eq. (12. 9 p reduces to Eq. (12. 3p . 

Next, we express the first two moments of the random variable Yt through the waiting 
time probability density p{t). Since the moments of N{t) are known from the CTRW 



theory (see, e.g., Ref. [16|), we reproduce here only the main results related to our situation. 
Introducing the probability P(n, t) that N{t) = n, we define the kth moment of the particle 
position Yf in the usual way: 

{Y,') = l'Y.n'nn,t) (2.10) 

n=l 

{k = 1,2,...). Then, using the Laplace transform of a function h{t), hg = C{h{t)} = 
dt e~^^h{t) (Res > 0), and taking into account that V{0,t) = drplr) and 

V{n,t) = f' dTp{T)V{n-l,t~T) (2.11) 
Jo 

{n > 1), we obtain 

{Yh. = l'^-^(ps^]'Y^- (2.12) 
s \ dps l-Ps 



Finally, applying to Eq. f l2.12p the inverse Laplace transform defined as h{t) = C ^{hs} 



f c+ioo 7 ^st 
'c-ioD 

we find the first 



{l/2ni) j^^l^ ds e'^^hs (c is chosen to be larger than the real parts of all singularities of /i^), 
:st 



and the second 



moments of Yt, which in turn determine the variance of the particle position: 

^\t) = {Y^) - {Y,)\ (2.15) 

III. CONDITIONS OF ANOMALOUS DIFFUSION 

As it follows from the waiting time probability density (12. 9p . the mth moment of the 
waiting time, = /g°° rfr r'^p(r) (m = 1, 2, . . .), can be written in the form 

7- = r i,ig'u(,)u(,') { ^^^^^if-^)" (3.1) 

J -go J -go \ f^[J +9) J 

Since the probability density u{g) is normalized on the interval [—go, go], from Eq. (13. ip it 
follows that < {ul + 2go/K,)"^/{f — go)"^- Thus, if f > go then all these moments are 
finite, and so in this case the classical central limit theorem for sums of a random number 
of random variables 



171] is applied to Yf. This implies that cr^(t) oc t as t — ?■ oo, i.e., 
at / > (70 the biased diffusion of particles is normal, and the rescaled probability density 
V{ip,t) = a{t)P{{Yt) + a{t)il),t) in the long-time limit tends to the probability density 
V{ip,oo) = (27r)~-'^/^e~'^^/^ of the standard normal distribution. 

It is clear from the above that the anomalous long-time behavior of the variance o''^{t) is 
expected at = oo when the mentioned central limit theorem becomes inapplicable. The 
condition = oo implies f = go that, in accordance with (12. 4p . yields Tmin = t^lf^go and 
Tmax = oo. Siucc the divergence of occurs when p(r) at r — )• oo tends to zero slowly 
enough, next we assume that 

P(^) - 7^ (3-2) 
(r —7- oo) with a > and a G (0,2]. Thus, the biased diffusion in a randomly layered 
medium is expected to be anomalous if both conditions, f = go and a G (0,2], hold. The 
former guarantees that the waiting time r can be arbitrarily large (rmax = oo), and so it is a 



necessary condition for anomalous diffusion. We note also that in this case one may expect 
that the rescaled probability density V{4j, t) at t — )■ oo approaches the stable probability 



density, as the generalized central limit theorem [18| suggests. 

The coefficient of proportionality a and the exponent a are in general not independent 
and can be found from the asymptotic behavior of the probability density u{g) in the vicinity 
of the point g = —go- Indeed, using the condition /f^^ dg gu{g) = 0, which is a consequence 
of the symmetry property of u{g), from Eq. (12. 9p at f = go we obtain 



p(r) ~ 4 



r2 



u\ — -go] (3.3) 



r 



uig)r^bl + ^\ , (3.4) 



(r — )• oo, K — )■ oo). Then, assuming that 

90 

where g — > —go, b > and /3 > 0, the asymptotic formula (13.31) takes the form 

BqA 1 



Pir)-bgo{-] (3.5) 



.90, 



Finally, comparing (13. 5 p with (13. 2p . we find a = /3 and 

a = aod(l-^^, (3.6) 

where Ood = ^fi'o(^Vfl'o)" is the parameter a in the overdamped limit. According to these 
results, the particle mass decreases the parameter a in comparison with the overdamped 
case but does not change the exponent a. 

An example of u{g) having the asymptotic behavior (13. 4p is the probability density 

r(/3 + l/2) g^V-' 

u{g) = U ~ ~ ' 3.7 

^ov^r(/3) V g^J 

where r(x) = Jl^ dy y^^^e~'^ is the gamma function, which corresponds to the symmetric 
beta distribution. According to Eq. (13.71) . this distribution is unimodal with the maximum 
at (yf = if /3 > 1, bimodal with infinite maxima at g = ±go if /3 < 1, and uniform if /3 = 1. 
Determining the parameter b directly from the density function (13. 7p . for this example we 
obtain 

r(« + i/2) f2uiY 

aod = ^ ^-p. s • (3.8) 

2V7rr(a) \ ^0 / 
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IV. LAWS OF ANOMALOUS DIFFUSION 



A. Long-time behavior of the inverse Laplace transform 

From a formal point of view, the moments (Yt) and (Y^) completely determine the vari- 
ance cr^(t). But the calculation of the inverse Laplace transforms in Eqs. fl2.13p and f l2.14p 
is a difficult technical problem because of the contour integration in the complex plane s. 
Fortunately, in the long-time limit this problem can be avoided. Such a possibility provides 



the celebrated Tauberian theorem for the Laplace transform [18| which is widely used in 
the theory of CTRW and its applications. According to this theorem, if h{t) is ultimately 
monotone and 



(0 < p < oo) as s — )■ then 



K = C{h{t)}^L[-^^^ (4.1) 



h{t) = C-\h,}^^^L{t)t^-' (4.2) 



as t — )■ oo, where L{t) is a slowly varying function at infinity. The term 'slowly varying' 
means that L{Xt) ~ L{t), i.e., lim^^oo -^('^^)/-^(^) = 1? for all A > 0. We note that, in 
contrast to Eqs. ( 12.13^ and f l2.14p . the parameter s in (14. ip is assumed to be a positive real 
number. 

The Tauberian theorem in the above form permits to find only the leading terms of the 
asymptotic expansion of the moments (Yt) and {Y^) as t — )• oo. If limt_j.oo(^t)^/(^t^) 7^ 1 
then, according to the definition fl2.15p . these terms completely determine also the leading 
term of the long-time expansion of cr^(t). However, if limt^^{Yt)'^ / (Y^) = 1 then for finding 
the leading term of cr^(t) at least the first two terms of each of the asymptotic expansions 
of (Yt) and (Y^) should be evaluated. Remarkably, these terms can also be determined 
from the Tauberian theorem if the Laplace transforms {Yt)s and {Y^^)s at s — )■ have the 
asymptotic form 

h.-^r^L^)!. (4.3) 
s'l \sJsP ^ ' 

with T] > p. In this case, replacing hg by hg — g/s^ and using the exact result C^^{l/s^} = 



t^~^/Y{ri) (loj, from (|42]) we obtain 
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It should be noted that qt^~^/T{ri) and L{t)tP~^ /V{p) actually represent the first two 
terms of the long-time expansion of h{t) only if rj > p. In the opposite case, when L{l/s)/s'^ 
is the leading term of the asymptotic expansion of hg, this may not be true. The reason is 
that in this case the second term of the asymptotic expansion of C^^{L{l/s)/s^} at t — )■ cxo 
may not be negligible in comparison with qf^'^ /T{ri). To illustrate this fact, let us assume 
that L{l/s) = k\n{l/s) {k is a scale factor) and p = 2, and consider the Laplace transform 
h, = q/s'' + (A;/s2)ln(l/s). Since £-^{{1/ s'^)\n{l/ s)} = t(lnt + 7 - 1), where 7 = 0.5772 



is the Euler constant 



19| . for the inverse Laplace transform of hg we obtain an exact result 



h{t) = qt^~^ /r{ri) + kt{lnt + 7 — 1). li rj > 2 and t — 00 then, keeping in h{t) the two 
leading terms, in accordance with (14. 4p we find h(t) — qf^^^ /T(ri) ~ ktlnt. However, if 
1] < 2 then the second term of the asymptotic expansion of k{'j — l)t, is 

not negligible compared to qt^~^/T{ri) and, as a consequence, the asymptotic formula (14. 4p 
does not hold. In this case only the leading term of h{t), kt In t, is determined from the 
Tauberian theorem. Thus, while at rj > p the first two terms of the asymptotic expansion of 
h{t) can be determined from the modified Tauberian theorem, Eqs. (14. 3 p and (14. 4p . to solve 
this problem in the opposite case it is necessary to go beyond the Tauberian theorem. 

Next, we use the Tauberian theorem, Eqs. (14. ip and (14. 2p . and its modified version, 
Eqs. (14. 3 p and (14. 4p . to find the long-time behavior of the first two moments, (It) and (Y^^), 
and the variance cr^(t) in the case of anomalous diffusion, i.e., when the conditions f = go 
and a G (0,2] hold simultaneously. Since the asymptotic solution of the CTRW is different 
for different intervals of a |20[, we consider the cases with a G (0, 1), a G (1, 2), a = 1, and 
a = 2 separately. 



B. ae{0,l) 

In this case it is convenient to represent the Laplace transform of the waiting time prob- 
ability density p(r) in the form 

/oo 
rfr(l-e~^Xr), (4.5) 
min 

which follows from the definition Ps = J^.^ dr e~*'^p(r) and the normalization condition 
/^.^^ dTp{T) = 1. Introducing the new variable of integration x = sr and using the asymp- 
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totic formula (13. 2p . we obtain 



l-Ps-s'^aj dx^-^ (4.6) 



as s — > 0. An integration by parts together with the integral representation of the gamma 
function |,21], r(a;) = dy e~'^y^~^, reduces (14.61) to the form 

aT(l — a) „ , , „, 

l-Ps ^— ^s". (4.7) 

Now, using this result and the Laplace transforms 

{Yt)s = I ^r^, {Y,')s = l'^^^, (4.8) 

s{i-Ps) s{i-psy 

of the first two moments of Yt, we find in the limit s — )■ 0: 

la 1 
ar(l - a) 



(yt)s - (4.9) 



and 

Since these asymptotic formulas are particular cases of the asymptotic formula (14. ip in which 
the slowly varying function L{l/s) is a constant, from (14. 2 p we obtain in the long-time limit 

<^->- .r(iJ;r(i + „) '° ("1) 

and 

Thus, in this case lim(_^oo(^)^/(^t^) 7^ 1 and the above asymptotic expressions yield 

^'^'^ ~ a-rkl- a) [nihaj ' POT^)'"" ^"''^^ 
According to this result, which agrees with that obtained in the context of the asymptotic 
solution of the CTRW [20], subdiffusion occurs if a G (0, 1/2) and super diffusion if a e 



(1/2, 1). If a = 1/2 then cr^(t) oc t and, in accordance with the commonly used terminology, 
the biased diffusion is normal. However, for normal diffusion processes both the mean and 
variance are proportional to time. Therefore, since (Yt) oc at a = 1/2, this type of 
diffusion should be more appropriately termed as quasi-normal. It is also worthy to note 
that, according to Eqs. (13. 6 p and (I4.13p . the larger is the particle mass, the stronger is 
diffusion. 
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C. ae{l,2) 



Since in this case Mmt^ooiYt)"^ / (y^) = 1 (see below), for finding the long-time behavior 
of cr^(t) we should determine the first two terms of the asymptotic expansion of (Yf) and 
(Y^) as t — )■ oo. To this end, taking into account that at a G (1, 2) the mean waiting time 
r = /^.^ dr rpij) exists, we use the following formula: 

("OO 

Ps = l-Ts- / dril -ST- e-'^)p{T). (4.14) 



Proceeding in the same way as before, we obtain 

aT{2 — a] „ , ^ 

Ps-l + rs } -fs" 4.15 

a{a — 1) 

as s — )■ 0, and the straightforward calculation of the Laplace transforms (14. 8 p yields 



/ 1 laVi2-a) 1 , 



and 



IY^\ 4/^«r(2-a) 1 

^ * r2 s3 ^ r3a(a-l) s^-"' ^ 

These asymptotic formulas are particular cases of the asymptotic formula (14. 3 P with 
L(l/s) = const. Therefore, using the well-known property of the gamma function, r(l+x) = 
xr(a;), from (14.41) we get 

(Yt)--t^^, ^^T7 7^'"" (4.18) 

and 

-^e^ ^- , t^-" (4.19) 

t2 r3a(a-l)(2-a)(3-a) ^ ' 

as t — )■ oo. Accordingly, the long-time behavior of the variance cr^(t) is described by the 
power law 

a\t)^^^, 77 ^t^'"- (4.20) 

Thus, since a G (1, 2), the transport of particles is super diffusive. Interestingly, depending 
on the exponent a, the increase of the particle mass fi can either enhance or suppress the 
biased diffusion. In order to show this, we first use Eq. (13. ip with m = 1 to represent the 
mean waiting time r in the form 

T = T„,(l + ^-^). (4,21) 
12 



where Tod = I~°godg u{g)/{go + g) is the mean waiting time in the overdamped hmit. 
Then, using this formula for r, the expression (13. 6p for the parameter a and the condition 
fi^Tmin ^ 1, we obtaiu 

^ = S^fl-^ + ^^^V (4.22) 

According to this result, the biased diffusion is enhanced by the particle mass at (3 — a;)rod > 
Sul/go and is suppressed at (3 — a)rod < Sul/go. In particular, if the probability density 
u{g) is given by Eq. (13. 7p then 

ul a -1/2 

Tod = — , (4.23) 

go tt - 1 

and so the former case occurs at a G (1, 3/2) and the latter at a G (3/2, 2). 



D. a = 1 

Here our starting point is the Laplace transform of the waiting time probability density 
p{t) represented as 



Ps 



roo / n \ 

l-aqs- dT{l - (p{r) - - , (4.24) 



where Qs = /^.^ dr {1 — e *^)/r^. The advantage of this representation is that the term aqs 
accounts for the asymptotic behavior of p{t), p(r ) ~ a/r^ as r —t- oo, in an explicit form. 



With the definition of the exponential integral |2l|, Ei{x) = dye ^/y, Qs can be written 



as 

qs = s^—^ + sE.iO (4.25) 

(^ = STmin). Since p(r) — a/r'^ = o(l/r^) (r oo), the integral term in Eq. (14.241) at s — )■ 

can be neglected compared to aos- Therefore, taking into account the asymptotic formula 
I I 

Ei{^) ~ ln(l/s) [2i], we obtain 

~ as In i. (4.26) 

Using this result and Eq. (14.81) for calculating the leading terms of the Laplace transforms 
at s — 7- 0, 

from the Tauberian theorem, Eqs. (14.11) and (14. 2p . we find the long-time behavior of the first 
two moments 

(>^*)--T^, {yt')--,X-r (4.28) 
a in I In t 
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It should be noted that similar asymptotic formulas for (Yt) and (Y^) were obtained in 



Ref. [20j. But because of the use of the waiting time probability density p(r) of a particular 
form, the asymptotic formulas derived in that paper do not depend on the parameter a. At 
the same time, as it was shown above, the parameter a contains an important information 
about the role of quenched disorder and particle mass. In particular, Eq. (13. 6p shows that 
the moments (14.281) increase with the particle mass. 

Since Imit^ooiXt)'^ / (Xt) = 1; ^oi finding cr^it) as t — )• oo we need to know at least 
the two leading terms of the long-time expansion of (Yf) and (Y^). In principle, using 
Eq. (I4.25P and the integral term in Eq. (I4.24p . we could easily find the asymptotic behavior 
of 1 — — asln(l/s) as s — )■ and, in this way, obtain the second terms of the asymptotic 
expansion of {Yt)s and {Y^)s. However, in contrast to the previous case, a straightforward 
application of the Tauberian theorem to this case does not provide a precise determination 
of the second terms of the asymptotic expansion of (Yf) and (Y^). As it was argued in 
Sec. IIV Al in order to find these terms it is necessary to go beyond the Tauberian theorem. 



E. a 



In this case we use the following representation for the Laplace transform of p{t): 



Ps = I - Ts - ars - 

where 



OO 1 

dr —(1 - ST - e~n 



min 



r3' 



2 



. - LE^iO. (4.30) 

This form of Ps explicitly accounts for both the finiteness of r and the asymptotic behavior 
of p{t). Since p{t) — a/r^ = o(l/r^) (r — i- oo), at s — > we can neglect the integral 
term in Eq. (14.29^ in comparison with ar^. This, together with the asymptotic formula 
r, ~ -(sV2) ln(l/s) (s ^ 0), yields 

- 1 + rs ~ -s^ln-. (4.31) 

2 s 

Using this result, the Laplace transforms of the first two moments of at s — )■ can be 
written as 

, , / 1 /a 1 , 1 , , 

iY.}.---,~^-in- (4.32) 
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and 

2 2/2 1 2fa 11 

~ ^-^In-. (4.33) 



Therefore, in accordance with the modified Tauberian theorem, Eqs. (14 ■3p and (I4.4p . we 
obtain 



/ . la 

and 



(Yt) - -t Int (4.34) 



(F/)-^t2~^tlnt (4.35) 

As a consequence, the long-time behavior of the variance is described by the asymptotic 
formula 

cx2(t) ^ Ljtlnt. (4.36) 
r 

The fact that a'^it) increases faster than t is in accordance with the asymptotic formula 
(I4.20p . Indeed, while approaches t, the coefficient of proportionality between (J^{t) and 
t tends to infinity as a -^^ 2 — 0. We note also that the ratio a/r^ at a = 2 is determined by 
the same Eq. (I4.22p . Therefore, if u{g) is given by Eq. (13. 7p then the larger is the particle 
mass, the weaker is diffusion. 

In conclusion of this subsection we would like to draw attention to the differences between 
our model and one-dimensional iterated maps which generate trajectories according to the 
rule Xn+i = Xn + F{x„). It is usually assumed {22! that F{x) is an antisymmetric, F{—x) = 
—F{x), and periodic, F{x + N) = F{x) {N is an integer), function. Due to these conditions, 
there is no drift, i.e., the quantity (x) = {xn+t — Xn) can be taken to be zero, where the 
angular brackets denote an average over a properly chosen set of initial conditions of x„ 
and t plays the role of the number of iterations. The same property, (x) = 0, holds also 



for maps perturbed by time dependent noise [23| and quenched disorder [2j] with zero 
means. Thus, these maps are unbiased and so the variance of Xn+t — Xn, i.e., (x^) — (x)^, 
reduces to (x^). Using the CTRW theory, in most cases it is possible to write this variance 
as the inverse Laplace transform and use the ordinary Tauberian theorem for finding its 



long-time behavior 22|, |23|]. In contrast, in our model (Yt) grows with time and although 
the moments (Yt) and (Y^) can also be represented as the inverse Laplace transform, the 
variance cr'^it) = (Y^) — {Yt)"^ cannot. Therefore, if limt-^oo (Yt)"^ / (Y^) = 1, the leading 
term of the asymptotic expansion of cr'^{t) as t — )■ 00 cannot be determined by applying 
the ordinary Tauberian theorem. We have solved this problem (for a 7^ 1) by using the 
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modified Tauberian theorem. It should also be noted that adding to F{x) a weak uniform 
bias breaks the sjTiimetry of the system and, as a consequence, leads to a time dependence 



of (x) ^5|. However, the biased maps considered in [25| do not exhibit anomalous diffusion 



at long times. 

F. Role of thermal fluctuations 

We complete our analysis with a qualitative discussion of the role of thermal fluctuations. 
These fluctuations can be accounted for by adding the thermal noise term to the right-hand 
side of Eq. f l2.ip . In this case some important conclusions can be drawn from the asymptotic 
behavior of the correlator K{x,y) = {[U{x) — f/(y)]^) as |x — y| — )■ oo, where U{x) is the 
random potential that corresponds to the random force g{x) = —dU{x)/dx. Since under 
thermal fluctuations particles can move in both directions, the random force and potential 
should be determined on the entire x-axis. Using the properties of g{x) and the continuity 
condition for U{x) with U{0) = 0, we obtain 

U{x) = -{x -nl)g^''^ + U{nl). (4.37) 

Here, x G [nl,nl + l), n = 0,±1,..., and U{nl) = -/ Es=o ^^'^ + ^^^"^ if n > and U{nl) = 
lY^^=_ig^^^ if n < —1. Taking into account that the random forces g^'^^ with different s are 
statistically independent, the definition of the correlator K{x,y) and Eq. fl4.37p lead to the 
asymptotic expression 

Kix,y)r^lal\x-y\ (4.38) 

(|x -y\-^ oo), where cr^ = dgg^u{g). 

The systems with K{x,y) oc \x — y\ have long been a subject of extensive study (see, e.g., 
Refs. jsl, [261 and references therein). A remarkable result obtained for these systems in the 
overdamped regime is that there always exists a threshold value ftr of the external force / 
in which the depinning transition occurs. This transition is characterized by vanishing the 
average particle velocity in the pinning state, when / < /tr, while in the depinning state, 
when / > /tr, particles move with a nonzero average velocity which strongly depends on /. 



According to [27(|, if g{x) is a bounded function then /tr < go (/tr — )• go as the temperature 
approaches zero) and the mean first-passage time in the depinning and pinning states is finite 
and infinite, respectively. Since the moments of the first-passage time can be associated with 
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the moments of the waiting time, we may expect therefore that at nonzero temperatures and 
f < go the exponent a and so the character of anomalous diffusion becomes depending on / 
(the diffusion behavior at f > go is expected to be normal). Specifically, with decreasing of 
/ from go to ftr the exponent a should also decrease from 2 to 1, and if / < ftr then a < 1. 
Of course, in order to find the dependence of a on / and the diffusion laws a quantitative 
consideration of the problem is needed. It is especially important because the effects arising 
from the joint action of quenched disorder and thermal fluctuations are often unexpected 
and even counterintuitive. 

V. CONCLUSIONS 

We have studied in the long-time limit the unidirectional transport of particles which 
occurs under a constant force in a randomly layered medium. The influence of the layers is 
modeled by a piecewise constant random force whose values in different layers are assumed 
to be independent and identically distributed with bounded support. We have reduced the 
problem of the unidirectional transport, initially formulated in the framework of the motion 
equation, to a continuous-time random walk on a semi-infinite chain. The main statistical 
characteristic of this approach, the waiting time probability density, is expressed through the 
probability density of the random force and particle characteristics, including the particle 
mass. By analyzing the dependence of the moments of the waiting time on the external force, 
we have formulated the conditions under which the biased diffusion exhibits the anomalous 
behavior. It has been shown that this behavior may occur only if the external force is equal 
to the boundary value of the random force. 

In order to find in the anomalous regime the long-time behavior of the first and second 
moments of the particle position, we have used the Tauberian theorem and its modified 
version allowing, in most cases, to determine the first two terms of the asymptotic expansion 
of these moments. Within this approach, we have found, with one exception, the explicit 
asymptotic formulas for the variance of the particle position, i.e., the laws of diffusion. The 
time dependence of the variance is completely controlled by the exponent describing the 
asymptotic behavior of the probability density of the random force in the vicinity of its 
minimum value. It has also been shown that, depending on the value of this exponent, 
the particle mass in the weakly underdamped regime can either enhance or suppress the 
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anomalous diffusion without changing its time dependence. 
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